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ABSTRACT 

We present the first comprehensive study of short-timescale chromospheric Ha variability in M dwarfs using 
the individual 15 min spectroscopic exposures for 52,392 objects from the Sloan Digital Sky Survey. Our 
sample contains about 10 3 - 10 4 objects per spectral type bin in the range M0-M9, with a total of about 206, 000 
spectra and a typical number of 3 exposures per object (ranging up to a maximum of 30 exposures). Using this 
extensive data set we find that about 16% of the sources exhibit Ha emission in at least one exposure, and 
of those about 45% exhibit Ha emission in all of the available exposures. As in previous studies of Ha 
activity (Lu a /L\, i) we find a rapid increase in the fraction of active objects from M0-M6. However, we find 
a subsequent decline in later spectral types that we attribute to our use of a spectral type dependent equivalent 
width threshold. Similarly, we find saturated activity at a level of L Ha /L}, i « 10~ 36 for spectral types M0-M5, 
followed by a decline to about 10~ 4 3 in the range M7-M9. Within the sample of objects with Ha emission, only 
26% are consistent with non-variable emission, independent of spectral type. The Ha variability, quantified in 
terms of the ratio of maximum to minimum Ha equivalent width {Rew), an d the ratio of the standard deviation 
to the mean (<7ew/(EW)), exhibits a rapid rise from M0 to M5, followed by a plateau and a possible decline 
in M9 objects. In particular, Rew increases from a median value of about 1.8 for M0-M3 to about 2.5 for 
M7-M9, and variability with Rew ^ 10 is only observed in objects later than M5. For the combined sample 
we find that the Rew values follow an exponential distribution with N(Rew) oc exp[-(7?Ew _ l)/2]; for M5-M9 
objects the characteristic scale is Rew~ 1 ~ 2.7, indicative of stronger variability. In addition, we find that 
objects with persistent Ha emission exhibit smaller values of Rew than those with intermittent Ha emission. 
Based on these results we conclude that Ha variability in M dwarfs on timescales of 15 min to 1 hr increases 
with later spectral type, and that the variability is larger for intermittent sources. Future studies using this large 
sample will address the variability timescales, the variability of other chromospheric emission lines (e.g., H/3, 
Ca II H&K), and the origin of the highest amplitude events. 

Subject headings: stars: magnetic fields — stars: flare — stars: late-type — stars: activity 



1. INTRODUCTION 

The study of magnetic activity in fully convective low mass 
stars and brown dwarfs (spectral types late-M, L, and T) has 
progressed in recent years from the question of whether these 
objects produce stable fields to a quantitative investigation 
of how the fields are produced and dissipated. Unlike the 
solar-type aVl dynamo, which operates in the transition region 
between the radiative and convective zones (the tachocline; 
lParker|[T955l) . objects later than spectral type M3 can only 
support a convective dynamo. Numerical simulations of such 
a mechanism are still at an early stage, but they suggest that 
large-scale axisymmetric fields can indeed be generated in 
fully convective objects, at least for conditi ons that roughl y 
correspond to mid-M dwarfs (M ~ 0.3 M Q : lBrowningf 2008 ). 
Thus, observational constraints on the scale, geometry, and 
dissipation of the fields are essential. 

Several observational techniques are now being used to 
address these questions, including Zeeman measurements 
in Stokes / and V, which probe the strength of the inte- 
grated surface fields and their large-scale topology, respec- 
tively (iReiners & Basrill2007l:lbonati et alJl2008UMorin et alJ 
2008), and activity indicators such as radio, X-ray, and Ha 
emission, which trace the dissipation of the field and hence its 
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strength and geom etry (e.g., IWest et all 12004 iBergerl 120061 
iBerger et alj|2009l) . The activity indicators also provide in- 
sight into magnetic heating, and their temporal variability can 
potentially trace the field properties on small scales that are in- 
accessible to the Zeeman measurements. The current results 
from Zeeman measurements point to a transition from mainly 
toroidal and non-axisymmetric fields in M0-M3 dwarfs to 
predominantly polo i dal axisymm etric fields in mid-M dwarfs 
dDonafi et all 120081: iMorin et alJ 120081) with field strength 
of ~ 0.1 - 3 1 kG (IReiners & Basrill2007t iDonati et alJ 120081: 
IMorin et"afll2008l) . Studies of activity indicators point to a 
rapid decline in X-ray and Ha activity (i.e., Lxfi a /L\, {) at 
about spectral type M7, and uniform radio lumi nosity (i.e., 
increasing L ral j/Lboi) at least to spectral type L3 (IBerger et alJ 
2009 and refer ences therein). The disparate trends may be 
due to a decoupling of the magnetic fields from the increas- 
ingly neutral atmospheres, or to a shift in the magnetic field 
configuration. 

While X-ray and radio observations provide powerful in- 
sight into the nature of the magnetic fields, Ha chromospheric 
emission, which traces gas at temperatures of ~ 10 4 K, is 
more easily accessible for large samples as a by-product of 
standard optical spectroscopic observations. In recent years, 
large spectroscopic samples of M dwarfs have become avail- 
able through dedicated studies and large-scale surveys such as 
the Sloan Digital Sky Survey (SDSS). These extensive sam- 
ples have led to several important results concerning chromo- 
spheric activity. First, the fraction of objects that exhibit Ha 
emission increases rapidly from ~ 5% in the K5-M3 dwarfs 
to a peak of ~ 80- 100% around spectral type M7, followed 
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by a subsequent decline to a few perc ent in the L dwarfs 
(Idizis et alJl200(ilWest et all2(Ml2008h . Second, while the 
level of activity increases with both rotation and youth in F-K 
stars, it reaches a saturated value of L-n a /L\, Q \ » 10~ 3,6 in M0- 
M6 dwarfs, followed by a rapid decline to Lu n /Lb n ] 1Q~ 5 
by spectral type LO dHawlev et all 119961: iGizis et all l2000t 
IWest et a02004l). and a brea kdown of the rotation-activity re- 
lation (iMartfn & Bouyl2002l) . Finally, a small fraction (< 5%) 
of late-M and L dwarfs have been serendipitously observed 
to exhibit Ha flares that reac h the saturated em ission levels 
found in the early-M dwarfs dLiebert et al.ll2003l) . 

To uniformly address the latter point — Ha variability — 
we recently carried out spectroscopic monitoring observa- 
tions of about 4 M4-M8 dwarfs with known Ha emission 
dLee et al.l l2009). With observations of about 1 hr per source 
and a time resolution of 5 - 10 min, we found that about 80% 
of these objects exhibit Ha variability on a wide range of 
timescales, ranging in amplitude from tens of percent up to 
a factor of about 5. Indeed, the timescale distribution for vari- 
ability "events" is nearly flat from 10 min to 1 hr, with fewer 
events on timescales below 10 min. The variability ampli- 
tudes follow an exponential distribution with a characteristic 
scale of Max(EW)/Min(EW)- 1 » 0.7. Finally, we found 
tentative evidence for increased variability with later spectral 

type- 
Here, we extend our study of Ha variability by three orders 
of magnitude using SDSS time-resolved spectroscopic data 
for a sample of 52,392 M0-M9 dwarfs (Knapp et al. 2009 
in preparation; hereafter Paper I). For the first time we take 
advantage of the individual 15 min exposures, with a typical 
number of 3 exposures per source. Using these data we can 
thus probe Ha variability on timescale of 15 min to ~ 1 hr. 
In this paper we focus on the Ha variability amplitudes and 
their relation to the Ha activity level (Lna/^boi) and spectral 
type. We also re-investigate the relation between Ha activity 
and spectral type using the individual 15 min spectra, rather 
than the SDSS pipeline-combined spect ra, which were used in 
previous work (IWest et al.ll2004l [2008b . In g2] we define our 
sample selection and describe the analysis method for measur- 
ing Ha equivalent widths. The results for Ha variability and 
activity are described in §|3] Future papers will focus on the 
variability timescales, the highest amplitude events, the vari- 
ability of other chromospheric emission lines (e.g., higher- 
order Balmer lines, Ca II H&K, He I lines), and the relation 
between chromospheric variability and various source prop- 
erties, such as age. 

2. SAMPLE SELECTION AND DATA ANALYSIS 

The details of the sample selection are described in Paper 
I, where preliminary characteristics of the sample are also de- 
scribed. Here we provide a summary of our selection tech- 
nique. We examined the combined spectra of all point sources 
in the SDSS DR7 database with z < 19.5, i-z> 0.2, and 
r—i > 0.5 mag (uncorrected for interstellar extinction). Us- 
ing a set o f spectral templa t es, inc luding the cool star tem- 
plates from lBochanski et alj (120071) . we determined the best- 
fit spectral type for each object and removed all objects with 
stellar spectra earlier than K, or with the spectra of galaxies 
or quasars. We further removed spectra determined to be de- 
ficient according to the following criteria: inadequate signal- 
to-noise ratio; poor sky subtraction; large wavelength regions 
with missing data; and contamination by light from a nearby 
object or by emission from diffuse interstellar ionized gas. 
The resulting list contains 59,427 stars. After removal of stars 



of spectral type K, L and T, and removal of stars with subd- 
warf spectra, we find a total of 52,392 M dwarf stars. 

Using the coordinates of these objects, and their plate and 
fiber identifications, we obtained the individual 15-min pro- 
cessed spectroscopic exposures. The total number of expo- 
sures is 205,823 with a range of 3-30 exposures per source. 
The distribution of the number of exposures is shown in Fig- 
ureQ] 

2.1. Equivalent Width Measurements 

To measure the Ha equivalent width (EW) in each individ- 
ual spectrum we follow a three-step process. First, we de- 
termine whether Ha emission is significantly detected. We 
define the minimum required Ha emission line significance 
as 3cr relative to the continuum level root-mean-square (rms) 
noise. This condition ensures that the Ha emission line is sig- 
nificantly detected regardless of the signal-to-noise ratio, and 
hence the apparent brightness of the source (a combination of 
its spectral type and distance). In Figure[2]we plot the median 
3cr threshold as a function of spectral type. For spectral types 
M0-M4 the threshold is about 1 A, but it then climbs rapidly 
to about 3 A for M6 and about 6 A for M9. This spectral type 
dependent threshold inevitably leads to different Ha detection 
statistics as a function of spectral type compared to previous 
analyses, which used a fixed threshold o f 1 A equivalent width 
from M0 to L0 (IWest et al.ll2004ll2008l) : see £3] We note that 
since the continuum flux in the Ha region declines with later 
spectral type, the use of a fixed equivalent width threshold al- 
lows weaker, and less significant "lines" to be counted as gen- 
uine detections. Our spectral type dependent threshold there- 
fore provides more robust statistics on the Ha properties at a 
fixed statistical significance. A total of 21,200 spectra exhibit 
Ha emission at > 3a significance (10.3% of all the available 
exposures). 

Second, we identify exposures with possible cosmic ray hits 
in the Ha line region using the standard SDSS pipeline flags. 
We find 1335 exposures with such flags, or about 0.65% of 
all spectra. From visual inspection, however, we find that in 
about half of the cases apparently genuine bright Ha emis- 
sion lines are flagged as cosmic rays. As a result, we perform 
a more robust cosmic ray rejection using one of the two fol- 
lowing criteria in addition to the pipeline flag: (i) the flux of 
a flagged pixel within ±2.5 A of the Ha line center is more 
than twice the maximum unflagged pixel flux, or (ii) more 
than half of the pixels in the Ha line and continuum region 
(6545 - 6595 A) are flagged by the pipeline. The former crite- 
rion allows us to distinguish cosmic rays with a narrow width 
from genuine Ha emission lines whose width is set by the 
spectral resolution. The latter criterion eliminates exposures 
with questionable flux measurements in the Ha line region. 
Using these additional criteria we reduce the number of re- 
jected exposures to 746 (0.36% of all exposures). 

To evaluate the statistical efficacy of this procedure in re- 
moving cosmic rays and in avoiding the removal of genuine 
Ha emission lines, we estimate the cosmic ray hit fraction in 
5 A windows located in regions of the spectra that are devoid 
of emission features. Using a random set of about 13,000 ob- 
jects that cover the full spectral type range of our study, we 
find that the expected cosmic ray hit fraction in the Ha line 
region is about 0.30%, in good agreement with the fraction 
determined by our cosmic ray rejection criteria. This indi- 
cates that our sample is not biased by cosmic rays, and fur- 
thermore that we are not missing a substantial fraction of the 
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large amplitude variability events due to cosmic ray rejection. 

Finally, to measure the Ha line equivalent width (EW) we 
simultaneously fit the Ha line with a Gaussian profile and 
the continuum region with a third-order Legendre polynomial. 
The uncertainty in the equivalent width (<tew) is determined 
over the same spectral region using the error spectrum asso- 
ciated with each exposure. For spectra that do not pass the 
3cr threshold, we use this value as an upper limit. These up- 
per limits are crucial for assessing the range of activity for M 
dwarfs that exhibit a mix of active and inactive spectra. 

Example spectra with no Ha emission, constant Ha emis- 
sion, and variable Ha emission are shown in Figure [3] 

3. RESULTS FOR Ha ACTIVITY AND VARIABILITY 
3.1. Detection Statistics 

In our sample of 52,392 M dwarfs we find 8,223 objects 
(15.7%) with at least one spectroscopic exposure that exhibits 
Ha emission. This is verified by the independent analysis 
presented in Paper I, which used the combined spectra. The 
detection fraction as a function of spectral type is shown in 
Figure [4] For comparison, we also show the detection frac- 
tions from lWesT et al. (2004) who used the pipeline-combined 
SDSS spectra and a fixed 1 A detection threshold. As in pre- 
vious studies, we find a rapid rise in the Ha detection fraction 
from a level of a few percent in the range M0-M3 to about 
20% in M4. However, beyond spectral type M4 our results 
diverge from previous studies. Namely, we find that the de- 
tection fraction rises to a peak of about 40% in M5-M6, and 
subsequently declines steadily to about 10% in M9. This is 
in contrast to the continued rise fro m M4 to M8 (with a peak 
fraction of nearly 100%) found bv lWest et all J2008). This 
divergence is largely due to our spectral type dependent de- 
tection threshold (Figure [2j. We demonstrate that this is the 
case b y applying our t hreshold to the distributions of Lu a /L\, \ 
from lWest et alj (120041) . using their equivalent width conver- 
sion factors. We find that the resulting reduction in active 
fraction is about 15-30%, and it explains the difference for 
spectral types M5-M7. We still find a difference of about 20% 
in the M8 and M9 b i ns. It is p ossible that the co ntinued rise 
found bv lWest et ai1(l2004l) and lWest etail ([2008) was due to 
further contamination caused by their use of a fixed equivalent 
width threshold. 

Of the 8,223 objects with at least one active spectrum, 45% 
exhibit Ha emission in all of the available exposures. The 
fraction of active objects as a function of the Ha detection 
fraction is shown in Figure This distribution does not take 
into account the fraction of objects in the sample with various 
numbers of exposures (Figure [TJ. For example, the peaks at 
Ha detection fractions of 1/3, 2/3, and 3/3 reflect the large 
number of objects with 3 exposures. Conversely, the low frac- 
tion of objects with Ha detection fractions of 80-95% simply 
reflects the small number of objects with a sufficiently large 
number of exposures to populate these bins (i.e., more than 5 
exposures). 

Even without correcting for this effect we find that for the 
objects with a small number of exposures (3-5), which oc- 
cupy well defined bins, there is a systematic decline in the 
fraction of objects as a function of increased Ha detection 
fraction. For example, the bins corresponding to 1 /3 and 2/3 
Ha detection fractions (dominated by objects with 3 expo- 
sures) exhibit a decline in the fraction of active objects from 
3.0% to 1.7%. Similarly, there is a systematic decline in the 
fraction of stars with Ha detections fractions of 1/4,2/ 4, and 



3/4 (with values of 1.1%, 0.5%, and 0.5%, respectively), and 
a decline for objects with detection fractions of 1/5 to 4/5 
from 0.7% to 0.2%. We note that the rise at 100% detec- 
tion fraction partly reflects the combined contribution from 
sources with a wide range of exposure numbers, but also ap- 
pears to be a genuine effect (some of the objects in this bin 
have steady Ha emission; see Figure 0). The decrease in 
the fraction of objects as a function of increased Ha detec- 
tion fraction may reflect a typical timescale for the variability, 
such that variability ona~ 15 min timescale is more common 
than on ~ 0.5 - 1 hr timescale. 

To appropriately normalize the distribution for the fraction 
of objects with different numbers of exposures, we divide 
each Ha detection fraction bin by the number of sources that 
can potentially populate this bin. For example, the bin at an 
Ha detection fraction of 0.5 is normalized by the total num- 
ber of sources with an even number of exposures, while the 
bin centered at 0.35 is normalized by the number of objects 
with exposure numbers divisible by 3. The resulting normal- 
ized distribution is shown in Figure [6] We find a relatively 
uniform distribution below an Ha detection fraction of about 
1/3 at a level of about 4%, followed by a decline to about 
1-2% for detection fractions of 0.4-0.95; the number of ob- 
jects with an Ha detection fraction of 1 remains unchanged 
by definition. 

These detection statistics alone already point to a significant 
level of variability in M dwarfs. Even if we make the most 
conservative assumption that all of the objects with a 100% 
Ha detection fraction are non-variable, the fraction of all M 
dwarfs with variable Ha emission is at least 8.6%. 

3.2. Activity 

Before we turn to a discussion of the Ha variability, we 
investigate the Ha activity trends. A comprehensive anal- 
ysis of M dwarf Ha activity as a function of spectral type 
based on SPSS p ipeli ne-combined spec tra was presented by 
IWest et all (12004 and lWest et all (|2008|) . In Figure [7] we plot 
Ln a /Lb i as a function of spectral type for each object. We use 
both the median and the maximum values, and find similar 
trends. The median Ln a /Lb \ value for spectral types M0-M4 
is roughly constant at 10~ 3 6 , followed by a decline to a value 
of about 10 ~ 4 3 in M7-M9 . These values are similar to those 
found bv lWest et al.1 (120041) . For the distribution of maximum 
Lua/Lboi we find a nearly constant value 10~ 3 5 for M0-M4, 
followed by a decline to 10~ 4 1 in M7-M9. Thus, the max- 
imum Lu a /Lb \ values follow the same trend as the median 
values, but the difference between the median and maximum 
values appears to increase with later spectral type, indicative 
of an increase in the Ha variability. 

In Figure[8]we plot the distribution of the median Lu a /L\, i 
values in separate spectral type bins. The shift in the peak 
of the distribution to lower values as a function of later spec- 
tral type is appar ent, and follows the same trend found by 
IWest et all d2004l) . However, unlike these authors we find a 
wider dispersion in the early spectral types, primarily as a re- 
sult of the significantly larger number of objects in our sam- 
ple. The standard deviation values range from about 0.2 to 
0.28 with no clear trend. 

3.3. Variability 

We now turn to the primary focus of our study — Ha vari- 
ability. We quantify the variability of each object using the 
ratio of maximum to minimum Ha equivalent width, Rew = 
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max(EW)/min(EW), as well as the ratio of the standard devi- 
ation to the mean equivalent width, Sew = 0ew/ (EW). In the 
case of /?ew we calculate a lower limit on the variability when 
at least one equivalent width upper limit is available, while 
for i7ew we tevat the upper limits as 3a detections; the result- 
ing values are thus lower limits on the variability. Within the 
sample of objects with Ha emission in at least one spectrum 
(8,223 objects) we find that about 26% are consistent with 
steady Ha emission within the uncertainties on the individual 
Ha equivalent width measurements. This fraction is relatively 
uniform across the full spectral type range of M0-M9, with a 
possible mild decline in the later spectral types. 

The distribution of /?ew values as a function of spectral type 
is shown in Figure [9] For each spectral type bin we sepa- 
rate objects with a 100% detection fraction from those with 
at least one upper limit. In both cases we find a clear rise in 
Rew as a function of spectral type, from a maximum value of 
about 5 for M0-M3 to about 10- 17 for M5-M8, followed 
by a possible decline in spectral type M9. Indeed, the objects 
that exhibit an Ha variability of > 10 are nearly all M5-M9 
dwarfs. For the objects with a 100% detection fraction we 
find that 1.2% have R EW > 5, while 0.1% have R EW > 10. 
For the objects with partial detections the fractions are 13.2% 
and 0.5%, respectively. These fractions, and the overall wider 
distribution of R E w values for objects with partial detections, 
indicate that objects with intermittent Ha emission produce 
higher variability ratios than those with persistent Ha emis- 
sion. To investigate this point in more detail, we plot the dis- 
tributions of /?ew values binned by detection fraction for all 
objects with 3 and 4 exposures (Figures [TOl and fTTTi. In both 
cases we find that the distributions for objects with partial de- 
tections are significantly broader than for objects with 3 /3 and 
4/4 detections. 

We further find that the R EW values follows an exponential 
distribution (Figure|9]). For the objects with a 100% detection 
fraction we find N(R E w) oc exp[-(R E y/- l)/0.7] for /?ew < 
4.5. This result is identical to what we previously found for a 
sample of abou t 40 objects, w hich r eached a m aximum value 
of R EW re 5 dLee et alj|2009b . The lLeeetal] d2009t) sample 
was too small to provide constraints on larger variability ra- 
tios. Here we find that for higher ratios C^ew ^ 4.5), the 
distribution appears to flatten to an exponential profile with 
N(R E w) oc exp[-(/?EW~ l)/2]. The distribution of /?ew lower 
limits is generally flatter, and appears to follow a single expo- 
nential profile with N{R E w) oc exp[-(7?E\v _ l)/2]. Since these 
are lower limits, we conclude that the overall distribution of 
R E w values for M dwarfs is an exponential with a characteris- 
tic scale of 7?ew ^ 2. 

The distribution of /?ew values for each spectral type is 
shown in Figure [12] The more pronounced tail of large /?ew 
values for spectral types > M4 is evident. In addition, we 
also find a systematic increase in the median of the distribu- 
tion, from about 1.8 for M0-M4 to about 2.5 for M6-M9. In 
Figure [13] we show the same distributions, but combined for 
spectral types M0-M4 and M5-M9. In both cases the distri- 
butions are well fit by an exponential profile with a character- 
istic scale of R EW - 1 pa 2 for M0-M4 and w 2.7 for M5-M9. 
This indicates that high amplitude variability is more common 
in late-M dwarfs. 

Finally, similar variability trends are found using (Tew. with 
a steady increase from a maximum value of about 0.4 for 
M0 to about 1 for M5-M7, followed by a decline to about 
0.7 for M8-M9 (Figure Q3). As in the case of /? E w, the 
distribution for partial detections is broader than for objects 



with 100% detections, which follows an exponential with 
N(a E w) oc exp(-CTEw/0.18). The distributions for each spec- 
tral sub-type are shown in Figure[l5]and exhibit an increase in 
(Jew from a median value of about 0.28 for M0-M3 to about 
0.4 for M6-M9. 

4. DISCUSSION AND CONCLUSIONS 

We study the variability of the Ha chromospheric emission 
line in M dwarfs using the individual 15 min spectroscopic 
exposures from SDSS for an unprecedentedly large sample of 
52,392 stars. The typical timescale probed by these data is 
15 min to 1 hr. About 16% of the sources exhibit Ha emis- 
sion in at least one exposure, and of those about 45% have 
Ha emission in all of the available exposures. Only 26% of 
all the objects with Ha emission are consistent with steady 
emission, spread relatively uniformly from M0 to M9. The 
detection fraction as a function of spectral type exhibits the 
known trend of a sharp increase from a few percent in M0-M3 
to tens of percent in later objects. However, unlike the results 
of previous studies based on SDSS pipeline-combined spec- 
tra and a fix ed equivalent width threshold of 1 A dWest et al.1 
12004 120081) . we find a peak detection fraction of ~ 40%, with 
a steady decline beyond M6. We attribute this trend to our use 
of a spectral type dependent detection threshold of ~ 3 - 6 A 
(3 a) in M5-M9 objects. 

The Ha variability, quantified in terms of R E w, exhibits a 
substantial increase with later spectral type. For M dwarfs 
as a whole, the distribution of /?ew values appears to follow 
an exponential with N(R E w) oc exp[-(/?EW _ l)/2]; the char- 
acteristic scale is about 2.7 for M5-M9. We also find that 
objects with partial detections exhibit a wider distribution of 
R E w values than those with persistent Ha emission. These 
results indicate that Ha variability increases with later spec- 
tral type, even as Ha activity declines, and that M dwarfs with 
intermittent Ha activity are more variable than those with per- 
sistent Ha emission. 

In the context of chromospheric heating, these trends sug- 
gest that the magnetic energy input typically varies by a factor 
of about 2 on timescales of ~ 15 min to ~ 1 hr. Moreover, the 
increased variability as a function of later spectral type, with 
a particularly large increase beyond spectral types M3-M4, 
hints at an overall shift in the magnetic field, leading to less 
uniform field dissipation. Taken at face value, this conclu- 
sion suggests that the relative contribution from small scales, 
which are likely to dissipate on short timescales, increases for 
later M dwarfs. This possible trend is in conflict with pre- 
liminary results from Zeeman measurements of a few objects 
in the range M0-M5, which point to an increase in the rel- 
ative contribution from large-scale fields dDonati et al.ll2008l: 
iMorin et al.ll2008h iReiners & Basrill2009h . Alternatively, the 
increase in Ha variability, particularly for sources with inter- 
mittent emission, may be due to an increase in the stochastic 
dissipation of the field in the presence of increasingly neutral 
atmospheres. In this scenario, heating of the bulk chromo- 
spheric plasma is suppressed by the decoupling of the field 
from the atmosphere (leading to a decline in Ln a /L\, {), but 
small-scale stochastic heating still takes place (leading to an 
increase in the Ha variability). The role of these two sce- 
narios may become clearer as samples of late-M dwarfs with 
Zeeman measurements become available. 
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FIG. 1. — Histogram of the number of M dwarfs in our sample as a function of the number of spectroscopic exposures. About 60% of the objects have 3 
exposures, 90% have 3-5 exposures, and 2.5% have > 10 exposures. 
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FIG. 2. — Median value of the equivalent width 3<r limit as a function of spectral type. The rise in the equivalent width threshold beyond M4 reflects the 
decrease in signal-to-noise ratio of the spectra due to a reduction in the continuum flux level. 
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FIG. 3. — Representative examples of M dwarfs with no detectable Ha emission in any of the individual exposures (left), constant Ha emission (center) and 
variable Ha emission (right). In each panel the time sequence evolves from darker to lighter shade of gray. 
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FIG. 4. — Fraction of M dwarfs with at least one exposure exhibiting detectable Ha emission as a function of spectral type. Our results are shown in black, 
while the results from West et al. 1 2004), based on SDSS pipeline-combined spectra, are shown in gray. The total number of objects in each spectral type bin from 
our sample is shown at the top. The divergence between the two distributions beyond spectral type M4 is primarily due to our use of a spectral type dependent 
equivalent width threshold (Figureff), as opposed to a uniform threshold of 1 A from M0-L0 in West et al. (2004); see £0] 
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FIG. 5. — Fraction of M dwarfs with detected Hct emission binned by the fraction of exposures in which Ha emission is detected. The peaks at 1/3, 2/3, and 1 
are indicative of the large number of objects with 3 exposures (FigureQ}- We do not show the large spike at zero detection fraction, corresponding to about 84% 
of our sample. 
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FIG. 6. — Same as Figure[5]but corrected appropriately for the number of M dwarfs in the sample that can contribute to each Ha detection fraction bin. 
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FIG. 7. — Ha activity level (LHa/i-bol) plotted as a function of spectral type. We use both the median (gray) and the maximum (black) Ha equivalent width for 
each object. The lines connect the median value for each spectral type bin, and the median values are listed at the top and bottom for the maximum and median 
£<Ha/A>ol> respectively. 
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FIG. 9. — Ratio of maximum to minimum equivalent width (Rew) for objects with detected Ha emission as a function of spectral type. Circles designate 
objects classified as active in every exposure, while triangles mark objects with at least one active exposure and using upper limits on the minimum equivalent 
width for the inactive spectra. The histograms for both samples are shown in the right panels. The light gray lines are an exponential fit to these samples. The 
distribution for 100% detection fraction appears to exhibit a break at Rew ~ 4.5. 
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FIG. 10. — Distribution of the ratio of maximum to minimum equivalent width (Rew) binned by Ha detection fraction for all objects with 3 exposures. The 
i?EW values for 1/3 and 2/3 detections are lower limits. The objects with partial detections exhibit higher amplitude variations. 
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FIG. 1 1 . — Same as Figure [Tol but for all objects with 4 exposures. 
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FIG. 12. — Distributions of the ratio of maximum to minimum equivalent width (Sew) for each spectral type. The median Rew value for each spectral type is 
listed on the right. The overall trend of increased variability with later spectral type is apparent. 
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FIG. 13. — Distributions of the ratio of maximum to minimum equivalent width (Rew) binned for spectral types M0-M4 (top) and M5-M9 (bottom). In both 
panels the gray line is an exponential fit with resulting characteristic scales of Rew - 1 ~ 2.0 (M0-M4) and rss 2.7 (M5-M9). 
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FIG. 14. — Ratio of the standard deviation to the mean equivalent width for objects with detected Ha emission as a function of spectral type. Circles designate 
objects classified as active in every exposure, while triangles mark objects with at least one active exposure and using upper limits on the equivalent width for 
the inactive spectra. The histograms for both samples are shown in the right panels. The light gray line is an exponential fit to the sample of objects with 100% 
detection fraction. 
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FIG. 15. — Distributions of the ratio of the standard deviation to the mean equivalent width for each spectral type. 



